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In QCD sum rules with external fields, the double dispersion relation 
^^ \ is often used to represent the correlation function. In this work, we point 

out that the double spectral density, when it is determined by successive 
CN ' applications of the Borel transformation, contains the spurious terms which 

^ ■ should be kept in the subtraction terms in the double dispersion relation. 

OO ■ They are zero under the Borel transformation but, if the dispersion integral 

^rr \ is restricted with QCD duality, they contribute to the continuum. For the 

f^ \ simple case with zero external momentum, it is shown that subtracting out 

2^ ' the spurious terms is equivalent to the QCD sum rules represented by the 

single dispersion relation. 
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The QCD sum rule |l[ is widely used in studying hadronic properties based on QCD 0. 
In this framework, a correlation function is introduced as a bridge between the hadronic and 
QCD representations. In the QCD side, the perturbative part and the power corrections are 
calculated in the deep space-like region (g^ = — oo) using the operator product expansion 
(OPE), which is then used to extract the hadronic parameter of concern by matching with 
the corresponding hadronic representation. 

In matching the two representations, it is crucial to represent the correlator using a 
dispersion relation. Usually in the nucleon mass sum rule as an example, the single- variable 
dispersion relation is used. With this, the QCD correlator calculated in the deep space- 
like region can be related to its imaginary part defined in the time-like region, which is 
then compared with the corresponding hadronic spectral density to extract the hadronic 
parameter of concern. The hadronic spectral density contains contributions from higher 
resonances as well as the pole from the low-lying resonance of concern. To subtract out 
the continuum, QCD duality is invoked above a certain threshold where the continuum 
contribution is equated to the perturbative part of QCD. This duality restricts the dispersion 
integral below the continuum threshold in the matching. Therefore, the predictive power of 
QCD sum rules relies heavily on the duality assumption. Indeed, in the quantum mechanical 
examples, the parton-hadron duality works well for two-point correlation functions [Q. 

Often, within the QCD sum rule framework, a correlation function with an external 
field is considered to calculate for examples pion-nucleon couplings 0-|^, nucleon magnetic 
moment 0. In such a case, as the two baryonic lines propagate through the correlator at 
the tree level, the double-variable dispersion relation P] can be invoked to represent the 
correlation function. Namely, 



oo 



n(Pi,P2) = / dsi I ds2-, 2\7 2\ + subtractions . (1) 

The subtraction terms serve to eliminate infinities coming from the integral. They are 
usually polynomials in pi or P2, which vanish under the Borel transformations. Thus, the 
subtraction terms should not contribute to the sum rules. As before, QCD duality is imposed 
to the correlator, which restricts again the dispersion integral below a certain threshold for 
both integration variables, si and S2- 

In general, the double spectral density p(si, S2) in the double dispersion relation of Eq. (|I]) 
is obtained formally by matching the correlation function with its corresponding OPE 11°^'^ [| 
calculated in the deep Euclidean region using QCD degrees of freedom. That is, 

II{pI, pD = n°P"(p2, pD + subtractions . (2) 

The LHS contains the spectral density in the integrals as given in Eq. (Q). To solve for 
the spectral density p{si, S2) from this formal equation, successive Borel transformations are 
needed to apply [^,0] on both side. This will eliminate the unnecessary subtraction terms. 
In doing so, the integrals are disappeared, reducing to a simple equation for the spectral 
density in terms of a given OPE. However, such a chosen spectral density, when it is put 



^In this work, we focus mainly on the OPE terms which contribute to the continuum. 



back to the double dispersion integral Eq. (|T]), can reproduce the original OPE up to some 
subtraction terms. Of course, the additional subtraction terms do not matter if a sum rule are 
constructed using precisely Eq. (0). But in fact, in constructing continuum contribution, the 
duality argument is imposed, which further restricts the dispersion integral. If this restricted 
form of the dispersion integral is used, the subtraction terms do not vanish even after the 
Borel transformations. In this work, we point out with some explicit examples that QCD 
sum rules using the double dispersion relation contain these spurious contributions. 

To proceed, we first demonstrate how the spectral density in the double dispersion rela- 
tion is usually determined 0|T^. The Borel transformation B{M'^,Q'^) is defined as 



Q^,n-*oo, Q2/n=Af2 nl \ dQ^ J 

With this definition, the Borel transformation converts the Q^ dependence of the function 
/ into the Borel mass dependence, M^. In doing so, polynomials in Q^ vanish. By applying 
the double Borel transformations on Eq. (|l|), we obtain, 

B{Ml-pl)BiMl-pi)Uiplpl) = dsj ds^pis^s,) e-^^IMl-s,iMl ^ ^^^ 

JO JO 

where we have used the formula. 

To eliminate the integral, we further perform additional double Borel transformations and 
obtain, 

B{tI -^)B{rl -^)B{Ml -pl)B{Ml -plmplpl) = p(-l, 1) , (6) 

where another formula for the Borel transformation, 

B{M^, Q2)e-'^'«' = M^6{a^M^ - 1) , (7) 

has been used. Note in this derivation, Si = l/rf > and S2 = 1/t| > 0, thus restricting 
the spectral density only in the region Si,S2 > 0. 

The OPE spectral density can be obtained by applying this operation on a given OPE. 
Note also that the integral interval should include the point provided by the delta functions. 
If the interval does not include the point, Si = 1/rf or S2 = 1/t|, then Eq. (|^) is not 
conclusive. We stress that the spectral density determined via Eq. (^) is correct within 
this context, obtained by successive applications of Borel transformation on the dispersion 
integral limited from zero to infinity. The subtraction terms, as they vanish under the Borel 
transformations, can be chosen freely as we like. 

In practice, however, QCD sum rules require a certain assumption for high energy part 
of the correlator, QCD duality. With this assumption, the dispersion integral is restricted 
below the continuum threshold S'o, and the Borel-transformed sum rule becomes, 

I'" ds, r ds^p'^^'^s,, S2) e-^^l^'^l-^-l^'^l = {'" ds, ['" ds2P^'"'^s^, S2) e-^^l^'^l-^-l^^l (8) 

JO JO JO JO 



pP'^™(si,S2) is obtained from hadronic representation of the correlator while p°^'^{si,S2) is 
obtained via Eq. (|]) for a given OPE. The LHS restricted below the continuum threshold 
5*0 can be calculated directly using the spectral density obtained from Eq.(^. Another 
equivalent method which is more useful for our discussion is to calculate the LHS via 



B{Ml-pl)B{Ml-pi) 



J Sn J Sn 



P°P^(S1,S2) 



So J So {S1-PJ){S2-P2) 



(9) 



Note that the continuum is subtracted out from the OPE using the duality argument. 
Integral intervals like J^ dsi /q ° ds2 or /g ° dsi J^ ds2 do not contribute because, as we will 
see, the spectral density is of the form ~ 5{si — S2) at least in our examples that will be 
considered in this work. The integral in the second term is bounded below by 5*0. As we 
have stressed above, in determining the spectral density via Eq. (^, it is important that 
the dispersion integral is limited from zero to infinity. Since the second integral is bounded 
below by Sq due to the duality argument, it is not clear if the subtraction terms as written 
in Eq. (|^) do not participate in the sum rule. This is our main question to be addressed in 
this work. 

Let us proceed how our question is realized in QCD sum rules with external fields. To 
do so, we consider as an example the two-point correlation function with a pion. 



n(g,p^) =ij d^xe^''-^(0|T[Jjv(x)Jjv(0)]|7r(p^)) , (10) 

where J^ is the nucleon interpolating field proposed by loffe O]. To be specific, let us take 



a typical OPE from this correlation function, 

■Qope _ / ^^ (pp{u) ln[—{q — uPt^Y] r^ I du (pp{u) ln[—upl — (1 — u)pf\ . (11) 



= / du fp{u) ln[—{q — uPt^^] ~ / du fp{u) In 
Jo Jo 


-upl 


q^ and p^ = {q — p-KY- We have taken the limit pi = 


--ml 



Here p\ = q and P2 = [q — P-k) ■ We have taken the limit p^. = ml = as is usually done in 
the light-cone QCD sum rules |]I2| . Note, Eq. (|TT]) contains terms polynomials in p\ or p\ but 
we did not specify these subtraction terms explicitly. For the twist-3 pion wave function, we 
take its asymptotic form, ^p{u) = 1 [|T^]. With higher conformal spin operators, the wave 
function takes more complicate form but our claims in this work are still valid even with 
more general wave functions. We will discuss this point later. 

To obtain the double spectral density, we take the operation as given in Eq. (^. For 
n°^*', we straightforwardly obtain 

pT^isi, S2) = -Si5{si - S2) . (12) 

Note, the spectral density is defined only in the region Si,S2 > 0. Therefore, the spectral 
density should be understood as being multiplied by the step function, 9{si)9{s2)- To 
include entire region of p{si,S2) 7^ 0, the lower boundary of the dispersive integral should 
be understood as 0~, infinitesimal but negative value. This subtlety does not matter in this 
example but is important in later examples. 

Normally, p'l"^ is simply used in QCD sum rules Eq. (|]) without justifying its use carefully. 
To see a problem with this spectral density, we put this expression into Eq. (P and perform 
the integrations using the Feynman parametrization. 
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S — UP2 — (1 — u)pI 



+ ln[s 



upl 



[s — UP2 — (1 — u)pi]'^ 

-il-u)pl] 



(13) 



The second term in the last hne is the anticipated logarithmic term matching the OPE of 
Eq. (^TJ). In other words, the second term is enough to reproduce the Borel-transformed 
OPE of Eq. (0). This means that the first term is spurious and it vanishes under Borel 
transformations with respect to the variables, —pf and —p^- Therefore, it is a part of 
subtraction terms and should not contribute to the QCD sum rule. That is, we have to 
subtract out this term using our freedom to choose any subtraction term. Of course, this 
subtraction by hand is not necessary if the sum rule is used in the context of Eq. (p. 
However, in practice, the sum rule is used in the context of Eq. (^ invoking QCD duality. 
The continuum part from this subtraction term. 
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(1 - u)pI 
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[1 - u)pI 



So 



becomes, under the double Borel transformation. 



.oMV-/- where ^^^ + ^. 



(14) 



(15) 



This nonzero continuum is spurious as it originated from the subtraction term. A prescrip- 
tion for including the continuum presented in Ref. [1^ can be obtained when this spurious 
continuum is kept and it is often used in the light-cone QCD sum rules |]T^. However, 
keeping this term while neglecting the OPE subtraction terms in Eq. (|lT]) is inconsistent. 
We now consider a slightly different OPE from Eq. (|TUp, 



UT 



duippiu) 



-upl 



u)pi 



(16) 



Again we take the asymptotic form for the pion wave function (pp{u) 



when the external momentum is zero, we have Hg^*^ = — 1/p^ [pf 



1 for simplicity. Note, 
pI = p"^)- It is clear 

that this OPE should not contribute to the continuum. Even if the sum rule is constructed 
with nonzero external momentum, this aspect should be recovered whenever we take the 
external momentum zero. Under the successive applications of the Borel transformation, it 
is straightforward to obtain the corresponding spectral density. 



pTisi,S2) 



S{si - S2)9{Si 



(17) 



Here we put the step function explicitly because the subtlety associated with the lower 
boundary affects the discussion. Using this spectral density, Eq. (||) becomes 



/ du 
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ds-, 



eis) 



upl~{l- u)plf 



(18) 



Note, we have for the lower limit of the integral in order to ensure that the integration 
includes entire region of p2(si, S2) 7^ 0. Integration by part leads to 



Jo- S — UP2 — (1 — U)Pi JO S — UP2 — (1 — U)pf 



(19) 



Here the first term yields the anticipated OPE of Eq. (0) and the second term, as the lower 
limit lies just below the zero, is zero. Note also that this separation becomes possible because 
the subtlety with the lower boundary. If there were no subtlety with the lower boundary, 
then we would not have the first term containing the delta function. If this were true, then 
in the limit of zero external momentum, Eq.(|l6D can not be equal to Eq . ([T9|) , which does 
not make sense. 

It is the second term that should be a part of the subtraction terms. Of course, this 
separation does not have any mathematical significance. It is however important physically 
because this separation enables us to identify where the spurious contribution to the sum 
rule comes from. That is, the second term, when the lower limit changes to 5*0, survives 
under the Borel transformations and contributes to the continuum. Once again, we have 
identified a spurious subtraction which contributes to the sum rule. 

Up to now, from the two simple examples, we have shown that QCD sum rules invoking 
the double dispersion relation contain the spurious terms originated from the subtraction 
terms. They contribute to the sum rule when QCD duality is imposed. This is a general 
statement as long as QCD sum rules are constructed using Eqs. (||),(|]) while the spectral 
density is determined via Eq. (|^). In general, for the correlator of Eq. (^) as an example, 
the OPE contains complicate functions like the pion wave functions. The general twist-3 
pion wave function can be written |T^, 



M^) = E ^ku'' . (20) 

k 

Using this general form, 112''*^ under the double Borel transformations becomes 

Muo)M^ (21) 

where uq = Ml/{Ml + M^) and M^ = M^M^/{Mf + M^). Under additional Borel trans- 
formations, the spectral density can be shown to be proportional to derivatives of 6{si — S2). 
The dispersion integral of Eq. (|T]) can be performed by integrations by part but, in doing 
so, the boundary terms become parts of the subtraction terms and produce the spurious 
continuum contributions when the sum rule is combined with QCD duality. In general, it is 
difficult to eliminate these spurious terms systematically. This is a generic problem in the 
sum rules using the double dispersion relation combined with QCD duality. 

Now, let us consider a simple case with the zero external momentum. As a specific 
example, we consider the two-point correlation function with a pion with the vanishing pion 
momentum (the soft-pion limit) or the same correlation function with one pion momentum 
taken out as an overall factor but in the rest with pt^ = |^,^,|l^ (beyond the soft-pion limit). 
Even in this case, the double dispersion relation is proposed as a correct representation of the 
correlator [Q. The double dispersion relation might be useful in treating phenomenological 
side properly but the spurious terms still persist. 

The double spectral density in this case takes the form 

p{si, S2) = p{si)5{si - S2) . (22) 



Since the two correlator momenta are equal in this case, the delta function appears as a part 
of the spectral density. The double dispersion relation Eq. (|l|) reduces to 



n(p2 
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subtractions 



Unlike to the single dispersion relation, the correlation function contains square of s 
the denominator. The spectral density p{s) is obtained by 



i3(r^^MM^V)^(p^) = ^ 



(23) 



p"^ in 



(24) 



=1/t2 



Thus, we can determine the derivative of the spectral density for a given OPE. The OPE 
corresponding to Eq. ( p^ is /n(— p^). We do not need to worry about the pion wave function 
(Pp{u) since its overall normalization, which is fixed to the unity, participates in the case. 



Substituting Ini^—p^) into Eq. (P^, we obtain 



dp{s) 



ds 



Pis) 



-s + constant 



(25) 



=1/t2 



The constant term, when put into the dispersion integral, yields the term 1/p^. To be 
consistent with the logarithmic behavior of the OPE, the constant should be zero. The rest 
of the spectral density leads to the dispersion integral. 



ds-— 
(s 



—s 



p" 



p 



,2^2 



S — p^ 



— ln{s — p^ 



(26) 



The second term in the RHS is what we have anticipated. This is what one would have 
obtained if the single dispersion relation is used. Again, this is enough to reproduce the 
Borel-transformed OPE of ln[—p'^). But, as before, the first term in the RHS is spuri- 
ous. This term is zero under the Borel transformation but the continuum gives nonzero 
contribution to the sum rule. Note that the first term can be separated as 



p^ 



pA 



p^ 



pz 



So 



+ 



p^ 



pz 



(27) 



So 



What is interesting is that the contribution from the upper limit in the first term cancels the 
one from the lower limit in the second term. These two terms coming from the continuum 
threshold survive separately under the Borel transform though their sum is still zero. If 
the spectral density of Eq. (|^) is simply used in the sum rule of Eq. (H), then the lower 
limit from the second term contributes to the continuum while the upper limit from the first 
term does not participate to the sum rule Eq. (H). Once again, the spurious nature of this 
continuum is obvious. 

Another way to support our claim is to consider the correlator Eq. (|1^) in the soft-pion 
limit. According to the soft-pion theorem, the correlator becomes the commutator with the 
axial charge Q^, 



n(g,p. = 0) 



rf^xe^«-^{ 



Q,,T[Mx)JNiO)] |0) 



(28) 



The commutator, if it is readily evaluated using the commutation relation for the quark 



fields, becomes the anticommutator. 



n(g,p. = 0) ~ |75, t J d^xe^^-^OlT [Jn{x)Jn{0)] |0)} . (29) 

It means, in the soft-pion limit, Eq. (|TUp is equivalent to the nucleon chiral-odd sum rule, 
which should be represented by the single dispersion relation in the construction of its sum 
rule. Note, in deriving this, we have used only the soft-pion theorem and the commutation 
relation for the quark fields. Therefore, this is an operator identity that has to be satisfied 
always. But if one starts from the double dispersion relation and takes the soft-pion limit 
afterward, then Eq. ( p9D is not satisfied exactly by the presence of the spurious term like 
in Eq. (|26|). 
Indeed, the spurious terms lead to different continuum factors as appeared in Ref. ||15|,|16 



P^ 



A similar discussion can be found in Ref. |T^ where we simply point out a pole at the 



continuum threshold, not clarifying its spurious nature. It is currently under-debate whether 



or not keep the terms in question in QCD sum rules with external fields [|T8|. But it is now 
clear from our discussion that they are spurious. Anyway, once the spurious subtractions are 
eliminated, then what is left is the same sum rule that we would have obtained via the single 
dispersion relation. Our argument can be generally applied to other OPE contributing to 
the continuum and it can be shown that subtracting the spurious continuum leads to the 
sum rule invoking the single dispersion relation. Normally, the continuum contribution is 
denoted by the factors, En{x = Sq/M^) = 1 — (l + a; + ----|- x"'/n\)e~^. Subtracting the 
spurious continuum replaces the continuum factor, -E„(x) -^ En-i{x) forn > 1 . The changes 
due to this spurious term are sometime huge as discussed in Ref. 0. In Ref. 0, figure 1 
shows for the sum rules with i'^^p^'^^ structure how this spurious continuum changes the 
Borel curve of irNN coupling. There, the Borel curve with the spurious continuum varies 
within the range 11.5 - 11.2. But without this spurious continuum, the variation scale 
becomes 15 - 30, clearly showing huge effects from the spurious terms. The ttNN coupling 
extracted from this Borel curve will be quite different from what we know experimentally. 
As discussed in Ref. 0, however, this only says that the Dirac structure i75P^7^ is not 
adequate for calculating the coupling. Instead, i'-y^ or "i^c^yq^p^ Dirac structure is more 
useful to calculate the coupling p. Nevertheless, from this consideration, we see that the 
spurious continuum sometime becomes substantial and thus should be treated carefully. 

In summary, we have pointed out in this work that QCD sum rules with external fields 
employing a double dispersion relation can contain the spurious continuum originated from 
the subtraction terms. This spurious term appears because the spectral density obtained 
via successive applications of the Borel transformation is not compatible with QCD duality. 
The spurious term should be subtracted out using the freedom for subtraction terms in QCD 
sum rules. In the case with the zero external momentum, subtracting the spurious term is 
equivalent to the sum rules invoking the single dispersion relation. Of course, our finding 
only affects continuum contributions but in some cases this modification leads to significant 
corrections to QCD sum rule results as is presented in Ref. 0. 
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